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Abstract 

Since it is not possible to predict when a Gamma-Ray Burst (GRB) will occur or 
when Active Galactic Nucleus (AGN) flaring activity starts, follow-up/monitoring 
ground telescopes must be located as uniformly as possible all over the world in 
order to collect data simultaneously with Fermi and Swift detections. However, 
there is a distinct gap in follow-up coverage of telescopes in the eastern U.S. region 
based on the operations of Swift. Motivated by this fact, we have constructed a 14" 
fully automated optical robotic telescope, Goddard Robotic Telescope (GRT), at 
the Goddard Geophysical and Astronomical Observator}'. The aims of our robotic 
telescope are 1) to follow-up Swift/ Fermi GRBs and 2) to perform the coordinated 
optical observations of Fermi Large Area Telescope (LAT) AGN. Our telescope 
system consists of off-the-shelf hardware. With the focal reducer, we are able to 
match the field of view of Swift- narrow instruments (20' x 20 ; ). We started scientific 
observations in mid-November 2008 and GRT has been fully remotely operated 
since August 2009. The 3o upper limit in a 30-second exposure in the R filter is 
~15.4 mag: however, we can reach to ~18 mag in a 600-second exposures. Due to 
the weather condition at the telescope site, our observing efficiency is 30-40% on 
average. 
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1 Introduction 


The 7 -ray emission which comes from extremely energetic 7 -ray sources such 
as Gamma-Ray Bursts (GRBs) and Active Galactic Nuclei (AGNs) (especially 
blazars) is believed to be produced by accelerated particles in a relativistic jet 
with bulk Lorentz factors of ~100 for GRBs and ~10 for AGNs. A GRB jet 
is formed when a massive star (> 10M©) collapses into a black hole (Woosley, 
1993; Paczyriski, 1998; MacFadyen & Woodley, 1999). On the other hand, 
an AGN jet is a continuous outflow from an active supermassive black hole 
(10 6-9 Mq) (e.g. Urry & Padovani, 1995). Although the energy scale and bulk 
motion of jets differ by an order of magnitude between GRBs and AGNs, 
their radiation processes are expected to be similar. Therefore, understanding 
radiation processes in the context of shock physics and particle acceleration 
using both GRBs and AGNs will provide a deeper understanding of the the 
fundamental physics in these extreme environments. 

Recent observations of prompt GRB optical emission by ground/space robotic 
telescopes are providing key data for understanding the radiation mechanism 
of GRBs. The observations of GRB 050820A by RAPid Telescope for Opti- 
cal Response (RAPTOR; (Vestrand et al., 2006)) and GRB 061121 by Swiff, 
UV/Optical Telescope (UVOT; (Page et ah, 2007)) indicate that there are 
at least two components in the prompt optical emission. One component is 
optical emission correlated with the prompt 7 -ray emission. The other com- 
ponent is a smoothly rising and decaying component during the prompt 7 -ray 
phase. The first component could be interpreted as emission from an internal 
shock because of the similar variability between the optical and 7 -ray bands. 
The second component could be due to an external shock interacting with the 
inter-stellar matter. However, the Robotic Optical Transient Search Experi- 
ment (ROTSE) has shown that there are a couple of cases where the early 
optical emission does' not correlate with the 7 -ray emission (e.g. (Rykoff et 
al., 2005)). The extremely bright prompt optical emission from GRB 0803 19B 
observed by TORTORA and ‘PI of the sky’ challenges the standard picture 
of the GRB emission model (Racusin et ah, 2008). 

Blazars form a sub-group of radio-loud AGNs and show an extreme vari- 
ability at all wavelengths (Urry, 1999). The most accepted scenario is that 
a rotating supermasive black hole surrounded by an accretion disk with an 
intense plasma jet closely aligned to the line of sight is responsible for the 
blazar emission. However, fundamental understanding of the radiation pro- 
cess in blazars requires extensive monitoring campaigns at all wavelengths. In 
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particular, properties of variability (including flares) and spectra from simul- 
taneous data in various wavelengths provides key information for a physical 
understanding of blazars. 


2 Motivation and Telescope Specification 


Since it is not possible to predict when a GRB will occur or when AGN flaring 
activity starts, follow-up/monitoring ground telescopes must be located as 
uniformly as possible all over the world in order to collect data simultaneously 
with Fermi and Swift detections. Based on the operations of Swift, however, we 
notice a distinct gap in follow-up coverage in the eastern U.S. region (Figure 1). 
This fact motivated us to construct a fully automated optical telescope at the 
Goddard Space Flight Center. If there is no GRB to observe, this telescope will 
perform extensive optical monitoring of Fermi Large Area Telescope (LAT) 
AGNs to identify the variability on the time scale of a day (Totani et al., 
2005). 

The Goddard Robotic Telescope (GRT) system consists of a 14" Celestron 
Optical Telescope Assembly (OTA), an Astro-Physics 1200GTO mount, an 
Apogee U47 CCD camera, a JMI electronic focuser (EV2CM) with a PC con- 
troller (PCFC), and a Finger Lake Instrumentation color filter wheel (CFW- 
1-8) with standard Johnson U, B, V, R, I and Clear filters. The quantum 
efficiency of the camera and the transmission curves of each filter are shown 
in Figure 2. With the focal reducer (Celestron f/6.3), a 20' x 20' field of view 
has been achieved. The observatory dome is an Astro Haven 7ft clam-shell 
dome. Opening and closing of the dome is controlled by the serial interface. 
Figures 3 and 4 are pictures of the telescope and the dome. 

The dome, the weather station (Davis Instruments Vantage Pro2), the we- 
bcams (Logitech Quick Cam Chat), and the rain sensors are controlled by 
another Linux PC. The dome is directly connected by a serial port to the PC. 
The Linux-based weather station control software (wview 1 ) and the webcam 
driver (spca5 2 ) are used to control the weather station and the webcams. 
The signal from rain sensors (Asuzac Inc.) is recorded by the I/O board (Ar- 
duino 3 ). and used to close the dome automatically in case of rain. 

The telescope site is the Goddard Geophysical and Astronomical Observatory 
(GGAO) which is about 1.5 miles northeast from the Goddard main campus. 


1 http://www.wviewweather.com / 

2 http: // mxhaard.free.fr/index.html 

3 http:// www.arduino.cc/ 
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Fig. 1. Distribution of the GRB follow-up telescopes (information from GCN Cir- 
cular of the Swift GRBs from September 2006 to March 2007). 
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Fig. 2. Left: The quantum efficiency of U47 CCD camera. The CCD chip is 
CCD47-10 (back illuminated CCD) by e2V technologies. Right: The transmission 
curves of each filter. 


3 Automation of the Observatory 


The whole telescope system, including the telescope mount, the focuser, the 
CCD camera and the filter wheel, is controlled by the ACP Observatory Con- 
trol Software . 4 An advantageous feature of the ACP software is to perform 
a short test exposure (five seconds in our case) after each slew and to do as- 

4 http://acp.dc3.com/index2.html 
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Fig. 3. Picture of the telescope. 


Fig. 4. Picture of the dome with 
the telescope inside. 


trometry automatically using the Guide Star Catalog 5 to find the center sky 
coordinate of an image. If the pointing is off, ACP requests the mount to re- 
slew so that the requested position will be exactly at the center of an image. 
Furthermore, ACP will automatically construct a pointing model by recording 
an offset between the requested position and the actual slew position. This 
pointing model will be used for every slew to achieve a better pointing at the 
initial attempt. 

We developed a script to make an observation schedule for a single night. 
This scheduler calculates the elevation of all objects in the database. It selects 
objects which are visible from the site (an elevation has to be larger than 35° 
to avoid the light pollution), and then schedules objects at the time when 
they reach the highest elevation. The output of the script is a text file with 
lines of ACP commands. ACP commands consist of an object name, a set of 
coordinates for each object, filters, exposure times of each filter, a bin size of 
images and the starting time of the observation. Since it is a text file, it is easy 
to modify the schedule by hand if we need to schedule an object which is not 
in the database. The sky flat, bias and dark frames are scheduled at the end 
of observations. The script runs once a day in the afternoon. This schedule 
file is automatically submitted through ACP by the observation control script 
which we describe next. As of September 2010, we have 58 Fermi LAT AGNs, 
70 variable stars, 33 near-by galaxies to search for supernova, and 37 guest 
observers’ objects in the database. 

A script to control the whole observatory has been developed. The script cal- 
culates the Sun's elevation every 10 minutes and sends commands to ACP 
or to the dome at specific Sun elevations to perform actions such as open- 
ing/closing the dome, changing the CCD temperature, submitting the schedule 
or re-submitting the schedule in case of a termination of on-going observation 
for various reasons. 


0 A user can specify either the Guide Star Catalog 1.1 or the USNO A2.0 catalog. 
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The GRB position information is received via a Gamma-ray bursts Coordi- 
nates Network (GCN) socket connection. Once the GRB position has been 
received, our custom GCN receiver code will calculate the elevation of the 
GRB at the site, and if it is observable (an elevation of a GRB has to be 
higher than 20° and also the Sums elevation has to be less than —20°), it will 
send commands to ACP to quit on-going observations and execute a GRB 
observation. We collect 200 frames of 5 sec, 80 frames of 30 sec, 120 frames of 
60 sec and 72 frames of 100 sec exposure of the GRB field in the time order 
given. Based on a real-time observation of several GRBs, GRT took 47 sec on- 
average to terminate on-going observations and complete a slew to the GRB 
location. The initial exposure was made 66 sec on average after we received 
the GRB position information through the GCN. 

The automatic script for LAT AGNs has also been developed to do the pho- 
tometry. The script performs the relative photometry using Sextractor 6 and 
PyRAF 7 adopting three reference stars defined for each object. The aperture 
photometry (IRAF phot task) is obtained using a circular extraction region 
of 5" for the source and an annulus with inner/outer radius of 5.5"/6.25" for 
the background. The script also extracts the publically available light curves 
of Fermi LAT, 8 Swift BAT, 9 and Swift XR,T 10 of the AGNs. Examples of 
multi-band AGN light curves produced by this script are shown in Figure 5. 
All of the AGN light curves observed by GRT are available to the public from 
the dedicated web page. 11 

In Figure 6, we summarize the automation of the observatory as a flowchart. 


4 Performance 


We used the observations of GRB 100915 to estimate the detection limit of 
our telescope system. The first several hours had good weather conditions. 
The CCD temperature was — 20°C. The observations were performed in the 
airmass range from 1.13 to 1.27 (about four hours) collecting 423 frames in 
total with four different exposures (5 sec, 30 sec, 60 sec and 100 sec). All the 
light frames are in the R filter and corrected for the bias, the dark and the flat. 
First, we extract the sources in the image using Sextractor. Then, ive perform 
photometry using the IRAF phot task. The aperture foreground radius is 1.6 


b http:/ /terapix. iap.fr/ rubrique.php?id_rubrique=91 / 

' http://www.stsci.edu/resources/software_hardware/pyraf 

8 http: // fermi.gsfc.nasa.gov/ssc/data/access/lat / msLlc/ 

9 http://swift.gsfc.nasa.gov/docs/swift/results/transients/ 

10 http://www.swift.psu.edu/monitoring/ 

11 http:/ / asd.gsfc.nasa.gov/Takanori.Sakamoto/GRT/index.html 
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Mrk 421 



Year 

Fig. 5. Multi-wavelength light curves of Mrk 421 (top) and 3C 279 (bottom). 
From top to bottom for each source, the light curve in the 0.1-300 GeV band by 
Fermi/ LAT, the 15-50 keV band by Swift/ BAT, the 0.5-10 keV band by Swift/XKT, 
and the optical R band by GRT. 


times the averaged FWHM of extracted sources. The background annulus 
with inner and outer radius are, respectively, 1.75 and 2.0 times the averaged 
FWHM of the extracted sources. A signal-to-noise ratio of extracted sources is 
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Fig. 6. Flowchart of GRT automation. 


calculated using the IRAF phot result. 12 The USNO B1 catalog R1 magnitude 
is used as the R magnitude of the extracted sources. The 3cr detection limit 
has been defined by fitting the signal-to-noise versus R magnitude data, and 
finding the 3cr crossing point to the fit (see Figure 7). The 3a detection limit 
of a 30 sec exposure in R is ~15.4 mag. To understand the detection limit as 
a function of exposure time, we repeated the same analysis of the images with 
different exposures (5 sec, 30 sec, 60 sec, 120 sec, 180 sec, 300 sec and 600 
sec). We stacked 30 sec images using swarp 13 to create images with >30 sec 
exposure. The relationship between the 3cr detection limit and exposure time 
is shown in Figure 8. We reached ~18 mag in a 600 sec exposure. 

In Figure 9, we summarize the fraction of good observations for each month 
for four AGN observations. Although we experienced few 7 months w r ith >60% 
of good weather night, the typical fraction of observable nights is 30-40%. 


12 The signal and the noise are calculated by 
signal = SUM - AREA x MSKY, 


and 


noise 


j SUM - AREA x MSKY 
V GAIN 


+ 1.23 AREA x ST DEV 2 . 


13 http://terapix.iap. fr/rubrique.php?id_rubrique=49 
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R; 30 s exposure 



Fig. 7. Relationship between signal-to-noise ratio and the R1 magnitude of the 
USNO B1 catalog for the image of the field of GRB 100915 at a 30 s exposure in 
the R filter. The horizontal dash-dotted line shows the 3cr detection limit. 
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Fig. 8. Relationship between 3cr detection limit and exposure time. The best fit 
relation (dashed line) is 3 cj upper limit = 14.1754 x log(exposure time)+13.9666. 

5 Current Status 



We started fully automatic operation in mid-November 2008. Although, at the 
beginning, we had conservative operations by having a duty operator every 
night, GRT is completely human free as of August 2009. We are currently 
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Fig. 9. Fraction of good observations for each month for four AGN observations: 
Mrk 501 (top left), Mrk 421 (top right), 3C273 (bottom left) and 1ES 1959+650 
(bottom right). 


observing ~85 objects and collecting ~350 images per night. 

Our optical contribution to the multi-wavelength observing campaign of the 
7 -ray outburst of 3C 454.3 has been published (Pacciani et ah, 2010). Further 
collaborations with the Fermi/LAT team for Mrk 421, Mrk 501, 1ES 1959+650 
and MG4 J2001 12+4352 are on-going. 

We have performed 13 follow-up observations of GRBs (12 Swift GRBs and 
one INTEGRAL GRB). Three of them were observed within 10 minutes after 
the GRB trigger. Although there is no detection of prompt optical emission or 
an afterglow so far, we confirmed that GRT can achieve a 17.5 to 19 magnitude 
limit in the R filter, and also respond automatically and promptly to a GRB 
through the GCN socket connection. 
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